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ABSTRACT 

We estimate the luminosity function of field galaxies over a range of ten magnitudes 
(-22 < M Bj < -12 for H = 100 km/s/Mpc) by counting the number of faint APM 
galaxies around Stromlo-APM redshift survey galaxies at known distance. The faint 
end of the luminosity function rises steeply at Mbj w — 15, implying that the space 
density of dwarf galaxies is at least two times larger than predicted by a Schechter 
function with flat faint-end slope. Such a high abundance of dwarf galaxies at low 
redshift can help explain the observed number counts and redshift distributions of faint 
galaxies without invoking exotic models for galaxy evolution. 

Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: 
luminosity function, mass function — cosmology: observations — surveys 

1. Introduction 

Knowledge of the space density of dwarf galaxies in the local Universe is of great importance 
both in its own right, for example to constrain models of galaxy formation and to understand 
the local distribution of matter, and also in order to interpret observations of faint galaxies. For 
quite some time now (see, for example, the review by Koo and Kron 1992), galaxy number counts, 
particularly in blue passbands, have been found to increase faster with apparent magnitude than 
predicted by simple no-evolution models, whereas the redshift distribution of galaxies in faint 
surveys is compatible with no evolution. Various models have been proposed to account for this 
discrepancy, including a disappearing or fading population of dwarf galaxies (Broadhurst et al. 
1988, Cowie et al. 1991, Babul and Rees 1992), non-conservation of galaxy numbers through 
merging (White 1989, Cowie et al. 1991, Broadhurst et al. 1992) or adoption of a non-zero 
cosmological constant (Fukugita et al. 1990). The Automated Plate Measurement (APM) galaxy 
number counts of Maddox et al. (1990c) show a factor of two higher counts at bj = 20.5 than 
expected from a no-evolution model, thus further exacerbating the discrepancy between observed 
number counts and redshift distributions. 
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Other authors, however, (Koo and Kron 1992, Koo et al. 1993, Gronwall and Koo 1995), have 
suggested that modifying the assumptions going into the no-evolution models can help reconcile 
the observed number counts and redshift distributions, without resorting to exotic evolution 
models or non-standard cosmologies. In particular, no-evolution models have traditionally assumed 
a Schechter (1976) form for the luminosity function with moderate faint-end slope, eg. a = —1.25 
(Ellis 1987) or a = —1.15 (Lilly 1993). Gronwall and Koo (1995) show that number counts in K, 
R and Bj bands, as well as colour and redshift distributions, are matched well by a no-evolution 
model in which the faint end of the galaxy luminosity function (Mb, ;> —150) is dominated by 
blue galaxies (B — V < 0.6) and rises significantly above a Schechter function with flat faint-end 
slope. To date, there have been few measurements of the shape of the local field galaxy luminosity 
function at such low luminosities. Eales (1993) measured an apparent upturn in the luminosity 
function (LF) at Mb, ^ —15 for galaxies at z < 0.1. However, he used the 1/V max method to 
compute the LF, and so it is subject to bias due to an inhomogeneous distribution of a very small 
number of galaxies. A similar caution should be applied to the results of Lonsdale and Chokshi 
(1993), whose lowest-luminosity measurement is anyway consistent with a flat faint-end Schechter 
function due to Poisson statistics alone. The analysis by Marzke et al. (1994) of the CfA Redshift 
Survey shows evidence for an upturn in the LF at low luminosities, although the possibility of a 
scale error in the Zwicky magnitude system makes the amplitude of the faint-end excess uncertain 
(but see Takamiya and Kron 1995). Probably the best evidence for an upturn in the LF at low 
luminosities comes from a recent measurement of the galaxy luminosity function from the ESO 
Slice Project (Zucca et al. 1997), who find an excess of galaxies fainter than Mb, ~ —17 above 
their best-fit Schechter function. 

In this paper we push the measurement of the field galaxy luminosity function to fainter limits 
than previous work by counting the statistical excess of faint (bj < 20.5) galaxies in the APM 
Galaxy Survey seen in projection around local galaxies at known distance in the Stromlo-APM 
Redshift Survey. The method by which we estimate the luminosity function is described in the 
following section. The galaxy samples and the application of the method are described in §H| and in 
§|] we justify the use of nearby centre galaxies to measure the faint-end of the galaxy LF. Results 
are presented in §[| and several tests of our procedure are described in §^. We conclude in §0. 

2. Method 

Given a catalogue containing positions and magnitudes for a large number of galaxies of 
unknown redshift and a smaller redshift survey in the same area of sky, one can measure the galaxy 
luminosity function to fainter luminosities than using the redshift survey alone. The technique 
used here to measure the abundance of dwarf galaxies relies on the assumption that correlated 
galaxies seen close in projection on the sky to a galaxy of known distance are also at the same 



1 Throughout, we assume a Hubble constant Ho of 100 km/s/Mpc. 
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distance. Although we do not know individually which galaxies are correlated, we can statistically 
determine the numbers of associated galaxies as a function of apparent magnitude, and hence 
absolute magnitude. Repeating this for a large number of centre galaxies, we can measure a 
luminosity function with good statistical accuracy at the faint end. 

Perhaps the most straightforward approach is to count the excess number of galaxies in bins 
of absolute magnitude out to a fixed projected separation around each centre galaxy, and to sum 
over centre galaxies. This approach was employed by Phillipps and Shanks (1987), who counted 
the numbers of galaxies down to bj = 20.5 from COSMOS scans of United Kingdom Schmidt 
Telescope (UKST) plates about centre galaxies taken from various Durham redshift surveys, to 
measure the field galaxy LF over the magnitude range —20.5 < Mbj < —16. However, as pointed 
out by Saunders et al. (1992) in a slightly different context, such an approach is not very efficient 
as centre galaxies at small distances contribute little signal but a large number of galaxies to the 
sum, since nearby a fixed projected separation corresponds to a large solid angle, and so the excess 
neighbours are dominated by projection effects. For this reason, Phillipps and Shanks only used 
centre galaxies at distances larger than 100/i _1 Mpc, thus limiting their ability to probe the faint 
end of the LF. 

Here, we generalize the method of Saunders et al. (1992) to estimate a (noisy) LF for each 
centre galaxy and then combine the LF estimates in a minimum-variance way. Saunders et 
al. counted galaxies in bins of projected separation a about centre galaxies of known redshift 
to estimate the projected correlation function, £(cr), which is related to the spatial correlation 
function, £(r), by 



We count galaxies as a function of a and magnitude to estimate the quantity X(M, a) = </>(M)H(<r). 
Assuming a model for H(cr) then yields the luminosity function 4>(M). 



Consider a centre galaxy at known distance y. We count the number of galaxies n in bins of 
projected separation a ± 5a/2 and apparent magnitude m ± 5m/2 about this centre galaxy. The 
expected value of n is 



where r 2 = x 2 + y 2 — 2xycos(6) is the square of the physical separation between galaxies at distance 
x and y with angular separation 9 = a/y, a is the fraction of the projected annulus within the 
survey boundary and 4>(M)5m is the number density of galaxies of absolute magnitude M ± Sm/2. 
The absolute magnitude M as a function of apparent magnitude m and distance x (in /i -1 Mpc, 
with corresponding redshift z) is given by the usual formula: M(m, x) = m — 5 lg[x(l + z)] — 25 — 3z, 




(1) 



2.1. Estimating X(M, a) from a single centre galaxy 





- 4 - 



where the final — 3z term is an approximate ^-correction for galaxies of unknown type in the bj 
passband. 

For centre galaxies at moderate distance y (y 3> a and <C 1), 

(n) rs N(m)5m + (f> \M(m,y) ]Smy 2 E(a) , (3) 

y L 

where N(m)5m is the surface density of galaxies of apparent magnitude m ± 5m/2 in the 2-d 
catalogue. 

As discussed by Saunders et al., this approximation is subject to biases due to the fact that 
for nearby galaxies, x 2 may grow faster than £(r) falls off, and so excess pairs may not be close to 
the galaxy of known redshift. Conversely, at large distances, the selection function may be falling 
off so steeply that neighbours in projection will on average tend to be nearer than y. These biases 
are corrected for as discussed below. 

To correct for the survey boundary, we count the number n r of randomly distributed points 
of mean surface density N r in the same a bins. Scaling n r to the surface density of galaxies of 
magnitude m, n'(m,a) = N ^ n r (a), we expect 



. . 27rasin0<5cr - 

(n ) = N(m)5m. (4) 

y 



Our estimate of the relative excess X(M, a) = ^(M)S(cj) is then 

1 



X(M, a) 



p(m,a,y) 

where 



n 



(5) 



p(m, a, y) = - k(m, a, y)y 2 (6) 
N(m) 

corrects for projection effects and the factor k(m, a, y) corrects for biases caused by the small 
angle/large distance approximation. From equation (||), the expectation value of X is 

1 f°° 

{X(M,a)) = ^- -/ (p[M(m,x)]x 2 £(r)dx (7) 



y 2 k(m,a,y) J 
and so setting 

S™<l>[M{m,x)WZ{r)dx 

k(m, a, y) = r- — (8 

y i =((T)0[M(m,2/)] 

will make X(M,a) an unbiased estimate of ^>(M)H(<t). Our estimator for 4>(M) thus depends 
on <j)(M). However if <j>(M) is modeled as a smooth function (eg. a Schechter function), a stable 
solution is reached within ten iterations. 
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2.2. Combining X(M, a) estimates to obtain <f>(M) 

Each Xi(M,a) calculated above for each centre galaxy should be an unbiased, albeit noisy, 
estimator of the quantity cf)(M)E(a). We now wish to combine these estimates in an optimal way 
to obtain a close to minimum variance estimate of <j>(M). 

Again following Saunders et al, the expected variance in the quantity rij is 

Var(n,) « nj[l + p(m, a, y)<p(M)E(a)} [1 + N(m)J 2 (9)][l + f<f>(M) J 3 (a)], (9) 

where 

J 2 (9) = 2vr / e'w(9')d9', (10) 
Jo 



and 



J 3 (a)=4vr f ' a'\{a')da' ', (11) 
Jo 



and / is the total fraction of galaxies with measured redshifts. In combining the individual 
measurements Xi to an overall estimate for </>(M), we sum over centre galaxies, weighting each 
measurement by bi = [H(<7)] 2 /Var(Aj), 

m) = aw^ (12) 

r l+p(m,a,y)cl ) (M)E(a)}[l + N(m)J 2 (9)}[l + f ( t } (M)J 3 (a)}. (13) 



6i n^[p(m,cr,2/)H((T)] 2 
The expected variance on our final measure of (f>(M) is then 



Var[0(M)] = !/]>>. (14) 



3. Galaxy Samples and Application 

For centre galaxies, we use the Stromlo-APM Redshift Survey (Loveday et al. 1996), which 
consists of 1787 galaxies with bj < 17.15 selected randomly at a rate of 1 in 20 from the APM 
Galaxy Survey (Maddox et al. 1990a,b). We count APM galaxies around each redshift survey 
galaxy in five linearly spaced bins of projected separation a up to 5/i _1 Mpc and in 22 evenly 
spaced bins in apparent magnitude from bj = 15 to 20.5. There are a total of 2,389,032 APM 
galaxies in this magnitude range. In order to correct for survey boundaries and holes drilled out 
around bright stars and large galaxies, we use a random catalogue of 3,476,477 random points 
uniformly distributed over the area of the APM survey. Although there are only about 1.5 times 
as many random points as galaxies overall, the number of random points is significantly larger 
than the number of galaxies (527,611) in the faintest apparent magnitude bin. Thus shot noise in 
the random catalogue introduces negligible random error into our results. 
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When counting galaxy and random points in projection around centre galaxies, we avoid 
objects within one major diameter of the centre galaxy. This is to avoid three biases that might 
otherwise be caused by a large foreground galaxy: (1) obscuration of faint galaxies directly behind 
it, (2) lensing of background galaxies and (3) break-up of a large galaxy by the APM machine 
into spurious sub-images, artificially boosting the apparent number of projected neighbours. The 
worst of these latter cases have already been removed by the "holes" drilled around large images 
during the construction of the APM survey (Loveday 1996, Maddox et al. 1990a). The majority 
of galaxies in the Stromlo-APM survey have jx^j ~ 25 isophotal major diameters of 20-100 arcsec. 

We assume a power-law form for the spatial correlation function, £(r) = (r/ro) -7 , with 
7 = 1.71 and ro = 5.1/i _1 Mpc (Loveday et al. 1995). The projected correlation function is then 
given by 

E(<r) = (a/a ) l -\ (15) 



where 

i 

7-1 



(16) 



For an initial luminosity function </>(M), we use the Schechter function fit from Loveday et 
al. (1992). We store the Xi calculated for each centre galaxy by projected separation bin and 
apparent magnitude bin in order to avoid incompleteness effects which would otherwise result if we 
used absolute magnitude bins. We convert apparent to absolute magnitudes when we combine the 
Xi estimates in equation (|T2|), using 22 bins in absolute magnitude from M = —22 to M = —11. 
At this stage, we also have the choice of which er-bins to use; for example we can count galaxies at 
projected separation 0-1, 0-2 or 2-5 /i _1 Mpc. The consistency of the estimates <j){M) for a range 
of projected separations will provide important confirmation of our results. 



4. Justifying the Use of Nearby Centre Galaxies 

Before we present the results of this analysis, we first justify our assumption that excess faint 
galaxies seen close in projection to nearby redshift survey galaxies really are physically correlated. 
The APM galaxy survey has a magnitude limit of bj = 20.5, corresponding to a limiting depth 
for L* galaxies of ~ 600/i _1 Mpc, whereas galaxies of absolute magnitude Mbj = —15 can only 
be seen out to a distance of 115/i _1 Mpc. One might therefore wonder whether fluctuations in 
the number of uncorrelated background galaxies might swamp the contribution from genuinely 
correlated dwarf galaxies for nearby centres. We therefore perform a simple experiment to show 
that the number of APM galaxies counted to various magnitude limits is systematically larger 
close to a nearby centre galaxy than in a random control field. 

Taking in turn those N cen = 477 galaxies in the Stromlo-APM Survey that are closer than 
100/i _1 Mpc as centres, we count the number of APM galaxies (iV ga i) within a projected separation 
of l/i _1 Mpc and also the number of background galaxies (Ab gr ) within the same solid angle at a 
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randomly chosen location. We use the random catalogue to correct for solid angle lost due to the 
survey boundary and holes drilled out around large images. This is done for APM galaxies to a 
magnitude limit of bj = 20.0, 19.0, 18.0, and 17.0. In figure |l|, we plot the frequency histograms 
(over the 477 centre galaxies) of the ratio N ga _i/Nb gr (solid histogram) and its inverse N^/Ng^ 
(dotted histogram). One can clearly see that the ratio Ng^/N^ is systematically greater than 
unity, ie. there is a statistical excess of APM galaxies near centre galaxies. One can quantify 
this excess by measuring the mean and rms of the ratio iVgai/^V"bgr over centres. The mean and 
its standard error [a ( iVg a i / )/\/ N ccn ] are shown in Table [p. Also shown in this table is the 
probability from the Kolmogorov test that the quantities N ga \ and N^ gr come from the same 
distribution. We see, particularly for the brighter magnitude cuts, that the distributions iVg a i and 
iVbgr are significantly different, and that the ratio N^/N^g,. is systematically larger than unity at 
at least the 4-cr level. From this test we conclude that we are justified in assuming that the excess 
APM galaxies seen close to Stromlo-APM centre galaxies are physically correlated, even for centre 
galaxies much closer than the characteristic depth of the APM galaxy survey. 

Table 1: Galaxy counts about centre galaxies iVg a i and control fields 



m fim 


(N gal /N hgr ) b 


K-S Test c 


a 20.0 


1.08 ±0.017 


2.3 x 10" 


-l 


b 19.0 


1.16 ±0.025 


4.9 x 10- 


-2 


c 18.0 


1.35 ±0.041 


2.2 x 10" 


-3 


d 17.0 


1.90 ±0.109 


2.1 x 10" 


-6 



a Magnitude limit of APM galaxies. 

b Mean and standard error in the ratio N ga \/Nb gI . 

c Probability that N ga ^ and N^ gI are drawn from the same distribution. 



5. Results 

In figure ||(a) the points with error bars show the first-iteration estimate of the luminosity 
function (p(M), counting galaxies to projected separation a < l/i _1 Mpc. For comparison, this 
plot also shows (as the continuous line) our earlier Schechter function fit to 0(M) using the 
redshift survey galaxies alone (Loveday et al. 1992) plus its extrapolation to lower luminosities 
(dotted line). Comparing with the earlier result, we see three regimes. (1) For M <^ —18, we 
see good agreement with the Schechter function. (2) For — 18 <^ M <^ —15 we see that the new 
determination of (f>(M) lies ~ 2a below the Schechter function. (3) At the faintest luminosities, 
M ;> —15, we see a sharp rise in (j>(M), rising well above the Schechter function. Thus there 
appears to be a significant excess of faint galaxies above a flat faint-end power-law, but we first 
need to understand why we disagree with our earlier estimate of (f>(M) over the magnitude range 
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-18£M£-15. 

We believe that this apparent discrepancy is due to the fact that we are really estimating the 
product 0(M)H(cr) and we are assuming that H(<r), and therefore is fixed, ie. independent of 
galaxy luminosity. In fact, as we showed in an earlier paper (Loveday et al. 1995), there is evidence 
that low-luminosity galaxies are less strongly clustered than L* (M* ~ —19.7) galaxies, and thus 
we would expect our estimate of 0(M)H(cr) to be biased low at low luminosities, as we indeed 
observe in figure ||(a) over the magnitude range — 18 <^ M <^ —15. This effect was also noted by 
Lorrimer et al. (1994) in their analysis of the distribution of satellite galaxies. The variation of 
clustering strength with luminosity is too poorly known to attempt to correct the <t>{M) estimate 
shown here. In fact, the present data provides probably our best constraints on how £(r) varies 
with luminosity. 

By absolute magnitude M ~ — 14 another effect is dominating: either the luminosity 
segregation reverses (such that M ^ —14 galaxies are more strongly clustered about L* galaxies 
than other L* galaxies) and/or the intrinsic space density of dwarf galaxies is significantly higher 
than predicted by extrapolation of a flat faint-end Schechter function. As we will see below, the 
latter seems a more likely explanation for the observed excess of faint galaxies near ~ L* galaxies. 
It also has the attractive feature of helping to explain the observed steep number counts of faint 
galaxies compared with traditional no-evolution models. 

Figure |||(a) showed our estimate of (f>(M) after just one iteration, ie. assuming that the true 
luminosity function is a flat faint-end Schechter function. Thus the steep faint-end is not due to an 
instability in our iteration procedure. Figure |||(b) shows our estimate of 4>(M) after ten iterations, 
by which time the solution has converged. We see that the faint-end slope has steepened slightly 
further. To model the observed 4>(M), we have fitted a modified form of the Schechter function, 
with an additional faint-end power law: 

«L, = ^(IO"e*p(^) 

In this formulation <p* , L* and a are the standard Schechter parameters, Lt is a transition 
luminosity between the two power-laws and (3 is the power-law slope of the very faint-end. No 
physical interpretation is intended by this choice of formula, it is merely a convenient way of 
modeling the observed 4>(L) over this extended range of luminosity and for estimating the faint-end 
slope. The line in figure |2](b) shows our best-fit "double power-law" luminosity function, which has 
parameters: a = -0.94, M* = -19.65, </>* = 1.54 x lCT 2 /i~ 3 Mpc 3 , M t = -14.07 and = -2.82. 
Although this fit is poor over the range — 17 <J, M <; —14, our 4>(M) estimate is almost certainly 
biased low over this range by the weaker clustering of galaxies fainter than L*. Clearly, the 
faint-end slope (3 = —2.82 cannot extend to indefinitely low luminosities, but it shows no obvious 
signs of flattening brightward of M = —12. 
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6. Checks and Tests 

6.1. Varying projected separation a 

In the preceding section we presented results based only on counts of neighbouring galaxies 
at a projected separation of a < l/i _1 Mpc. By increasing a one decreases shot-noise errors by 
increasing the total number of neighbours, although now a smaller fraction of these neighbours 
will be genuinely correlated with the centre galaxy. Figure || shows our estimate of 4>{M) if we 
count galaxies to a projected separation a < 2/i _1 Mpc (solid symbols) and a < 5/i _1 Mpc (open 
symbols). We see that the estimated 4>(M) is insensitive to the limiting a used, except that the 
deficit at M ~ —15 worsens as the maximum projected separation increases, and the error bars 
increase. This is to be expected from our earlier results on luminosity segregation, since figure 6 
of Loveday et al. (1995) shows that the weaker clustering of sub-L* galaxies is more pronounced 
at scales larger than l/i _1 Mpc. The error bars also increase if we decrease the limiting projected 
separation to values of less than l/i _1 Mpc. Counting neighbours to a ~ l/t — 1 Mpc, as shown in 
Figure g thus appears to be about optimal for this analysis. 

6.2. Varying limiting apparent magnitude 

The second test we have performed is to cut back on the magnitude limit to which we count 
APM galaxies. There are two reasons for doing this. The first is that there is some indication 
that star-galaxy separation is slightly less reliable in the faintest half-magnitude slice of the APM 
survey than for brighter objects (Maddox et al. 1996). The second is to show that the poor match 
between volumes sampled by the bright redshift survey and the faint photometric survey does 
not lead to a systematic bias in the estimated luminosity function. In Figure |] we plot the LF 
estimated by counting APM galaxies to 20th, 19th, 18th and 17th magnitude respectively. In 
each case, we see that the estimated LF is consistent with that estimated from the full bj < 20.5 
sample. Even cutting back the "faint" survey to bj < 17 (actually slightly shallower than the 
redshift survey), we still see evidence for a steep faint end to the luminosity function. One should 
note that the signal contributing to the faint-end of the LF is coming from decreasingly smaller 
volumes as one cuts back on the magnitude limit. A galaxy of absolute magnitude M# 7 = — 15 
can only be seen to distances of 93, 60, 39 and 25 /i -1 Mpc respectively for galaxy samples cut to 
bj = 20.0, 19.0, 18.0 and 17.0. The consistency of estimated faint-end slopes using galaxies over 
this range of volumes argues strongly against the observed steep faint-end slope being due to a 
"sampling" or "volume" effect. 
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6.3. Analysis of simulations 

As a further test of our methods, we have generated two sets of mock-APM surveys with 
known luminosity function and two-point correlation function, thus allowing us to check that we 
indeed are measuring the "true" luminosity function. The first set of simulations ("Sch") has 
galaxy luminosities drawn from a flat faint-end Schechter function (a = —1, M* = —19.5). The 
second set ("DP") has a "double power-law" function of the form ( |i~7| ) with parameters a = — 1, 
M* = -19.5, (5 = -2.6, M t = -14.5. The idea is to check that our method for estimating </>(M) 
can reliably distinguish between these two models. 

For each model we generated five Soneira and Peebles (1978) hierarchical clustering 
simulations, each containing 2.4 million galaxies in the area of the APM survey. From each of 
these mock APM surveys, we formed mock Stromlo-APM Redshift Surveys by sampling one 
galaxy in 40 brighter than bj = 17.15 at random. Note that the sparser sampling rate compared 
with the real data (1 in 20) is required to match the numbers of galaxies in the faint APM Galaxy 
Survey and the brighter redshift survey sample. This reflects the steep slope seen in the APM 
number counts by Maddox et al (1990c) but not modeled in the simulations. The mock redshift 
surveys contained on average 1648 galaxies each, and the measured real space galaxy correlation 
function was well described on scales 1-20 /i _1 Mpc by a power-law with parameters 7 = 1.8 and 
ro = 4.9 for the "double power-law" simulations, and with 7 = 1.8 and ro = 5.1 for the Schechter 
function simulations, thus providing an excellent match to the observed real-space clustering of 
Stromlo-APM galaxies (Loveday et al. 1995). 

We analysed the simulations in the same way as the real data; counting galaxies in the 
mock-APM survey about centre galaxies in the respective mock-Stromlo survey. In Figure ||(a) 
we plot the mean <fi(M) measured from the five "DP" simulations as the symbols; the error bars 
show the one-sigma scatter between the realizations. These error bars are in good agreement with 
the errors estimated from equation (|L1|). Also shown are the individual estimates (dashed lines) 
and the true luminosity function (continuous line). We see that the mean estimated <f>(M) is in 
reasonable agreement with the true LF: if anything we tend to underestimate the slope of the 
faint-end of the luminosity function. We also see that estimates from the individual simulations 
can differ significantly from the true LF: two of the estimates show a decrease in <f>(M) faintward 
of M ~ —17. This lack of robustness in the individual estimates of the LF is also apparent in 
Figure ||(b), where we plot the results from the "Sch" simulations. Although the mean 4>(M) is 
consistent with a flat faint-end, and inconsistent with the double power-law function shown, two of 
the estimates do show an apparent rise above the expected flat faint-end slope. Neither, however, 
maintains the rise as faint as M = —12. 

Overall, although individual simulations can show a wide deviation from the expected 
behaviour, it seems that we are more likely to underestimate rather than overestimate the 
faint-end slope of the luminosity function, and hence the space density of dwarf galaxies. Our 
estimated errors are in good agreement with the scatter between simulations. 
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7. Conclusions and Discussion 

We have seen that the number of dwarf (M ;> —15) galaxies seen close in projection on the 
sky to ~ L* galaxies (—22 <J M <^ —15) is much larger than expected for a flat faint-end Schechter 
function and for the standard galaxy correlation function (7 ~ 1.7, ro ~ 5/i _1 Mpc). We may thus 
infer that the space density of galaxies rises sharply above a flat faint-end Schechter function for 
M ^ — 15 and/or that the clustering of dwarf galaxies around ~ L* galaxies is much stronger than 
the auto-correlation function of ~ L* galaxies. An extreme case of the latter explanation might 
be that dwarf galaxies only exist close to ~ L* galaxies. 

It is no easy matter to determine which (or both) of these explanations is correct. Since all 
of the M ;> —15 dwarf galaxies must lie within y <^ 115/i _1 Mpc to be visible in the APM survey, 
even if there is a generally higher space density of dwarfs, they have a relatively minor (~ 20%) 
effect on total predicted number-counts in the survey. However, one can place a lower limit on the 
space density of dwarf galaxies by making the extreme assumption that they only exist close to 
~ L* galaxies. In fact, in Figure ||, we see an excess of dwarf galaxies up to a projected separation 
of at least 5/i _1 Mpc from ~ L* galaxies. Thus a lower limit may be estimated by assuming that 
they occur with the measured space density only within 5/i~ 1 Mpc of an ~ L* galaxy. The lower 
limit on the average mean density of dwarf galaxies is then given by multiplying their measured 
space density by the fraction of space within 5/i _1 Mpc of an ~ L* galaxy. In order to estimate 
the filling factor of ~ L* galaxies in the local Universe, we have generated another set of five 
Soneira-Peebles simulations within a cubic volume 100/i -1 Mpc on a side and without applying a 
selection function. These simulations contain 47,000 galaxies each, and thus their space density is 
n = 4.7 x 10 -2 /i 3 Mpc~ 3 , as measured for ~ L* galaxies in the Stromlo-APM survey (Loveday et 
al. 1992). Additionally, the two-point correlation function of galaxies in these simulations is set 
to be well fit by a power-law with 7 = 1.82 ± 0.06 and ro = 5.0 ± 0.2, in good agreement with 
the clustering measured by Loveday et al. (1995). Given these simulations, it is then simple to 
perform a Monte-Carlo calculation of the volume inside each cube within 5/i _1 Mpc of an ~ L* 
galaxy. We find that this factor lies in the range 0.57-0.62, ie. a randomly chosen point in space 
has a 60% chance of lying within 5/i _1 Mpc of one or more ~ L* galaxies. Now, integrating the 
luminosity function plotted in Figure ||(b) between M = —15 and —12 yields a measured density 
of dwarf galaxies of n ~ 0.20/i 3 Mpc~ 3 . Correcting this by the extreme assumption that dwarf 
galaxies are only found within 5/i _1 Mpc of an ~ L* galaxy results in a lower limit on the space 
density of dwarf galaxies of fi ~ 0.12/i 3 Mpc~ 3 . This is a factor of two higher than the density 
fi ~ 0.058/i 3 Mpc -3 inferred from the a = —1.11 Schechter function fit by Loveday et al. (1992). 

In fact, the true space density of dwarf galaxies is likely to be significantly higher than this, 
for a number of reasons. First, our estimator assumes that galaxy clustering is independent of 
luminosity, whereas we know that sub-L* galaxies are less strongly clustered than more luminous 
galaxies (eg. Loveday et al. 1995). If this luminosity segregation extends to dwarf galaxies, then 
we will have underestimated their space density. Second, analysis of simulations (§ |6.3D shows 
that our estimator tends to underestimate the faint-end slope of 4>(M) slightly, possibly due to a 
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Malmquist-type bias. Third, as discussed by numerous authors, most galaxy surveys are likely to 
be missing a substantial fraction of low surface brightness galaxies, many of which will be dwarfs. 
For example, Sprayberry et al. (1997) find a pronounced upturn in the luminosity function for 
their sample of low surface-brightness galaxies. Thus we regard our above estimate of the space 
density of dwarf galaxies, n ~ 0.12/i 3 Mpc -3 , as a lower limit on the true value. 

A high space density of dwarf galaxies, assuming that they are predominantly late-type, blue 
galaxies, which suffer smaller i^-correction dimming than redder, early-type galaxies, provides a 
natural explanation for the steep observed number counts of faint galaxies. Evidence for a large 
contribution from late-type galaxies to the faint galaxy counts in the HST Medium Deep Survey 
has been presented by Driver, Windhorst and Griffiths (1995). Zucca et al. (1997) have found 
that the luminosity function for emission-line galaxies (ELGs) is significantly steeper at the faint 
end than the LF for non-ELGs, also supporting the hypothesis that the faint-end of the galaxy 
LF is dominated by late- type galaxies. Gronwall and Koo (1995) were able to match observations 
of galaxy number counts in the K, R and Bj bands, as well as colour and redshift distributions, 
for a mild evolution model by assuming that the faint end of the galaxy luminosity function is 
dominated by blue galaxies (B — V < 0.6), and rises significantly above a Schechter function with 
flat faint-end slope. We plot the Gronwall & Koo model total luminosity function in Figure §(b), 
and see remarkably good agreement with our observations. Our results thus support the model of 
Gronwall & Koo; there is no need to invoke exotic forms of galaxy evolution to explain observed 
galaxy number counts at faint magnitudes. 

It is a pleasure to thank my colleagues George Efstathiou, Steve Maddox, Bruce Peterson 
and Will Sutherland, who made the APM and Stromlo-APM surveys possible. The simulations 
used in this paper were carried out using computing facilities at the Fermi National Accelerator 
Laboratory. The referee, Chris Impey, is thanked for his comments which have (hopefully) helped 
make these results more convincing. 
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Figure Captions 

Fig. 1. — Plots of the ratio of APM galaxies counted around centre galaxies (iV ga i) over the number 
of APM galaxies in randomly placed control fields (N^ gT ) of the same solid angle. The solid 
histogram shows the ratio Ag a i/iVbg r and the dotted histogram shows the inverse N^/N^. APM 
galaxies were counted to a bj magnitude limit of (a) 20.0, (b) 19.0, (c) 18.0 and (d) 17.0. 

Fig. 2. — Estimates of the galaxy luminosity function 4>(M). (a) After just one iteration. The 
continuous line shows the earlier Schechter function fit to <$>{M) using the redshift survey galaxies 
alone (Loveday et al. 1992); its extrapolation to lower luminosities is shown by the dotted line, (b) 
After ten iterations, by which time the solution has converged. The smooth curve shows a "double 
power-law" Schechter fit (eqn. |i~7| ) and the histogram shows the Gronwall & Koo (1995) model. 

Fig. 3. — Estimates of the galaxy luminosity function <$>{M) obtained by counting APM galaxies 
to a projected separation a < 2/i _1 Mpc (solid symbols) and a < 5/i _1 Mpc (open symbols). The 
smooth curve is the same "double power-law" Schechter fit shown in Figure |2|b. 

Fig. 4. — Estimates of the galaxy luminosity function <$>{M) obtained by counting APM galaxies 
to a bj magnitude limit of (a) 20.0, (b) 19.0, (c) 18.0 and (d) 17.0. The smooth curve is the same 
"double power-law" Schechter fit shown in Figure ||b. 

Fig. 5. — Analysis of simulations. The dashed lines show the individual realizations, the symbols 
with error bars show the mean and 1-sigma fluctuations between them. The continuous line shows 
the "double power-law" LF used in the "DP" simulations, (a) Results from the "double power-law" 
simulations ("DP"), (b) Results from Schechter function simulations ("Sch"). 
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